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ABSTRACT 

'<6: 

We investigate outflows driven by a giant protoplanet using three-dimensional 
MHD nested grid simulations. We consider a local region around the protoplanet 
in the protoplanetary disk, and calculate three models: unmagnetized disk model, 
magnetized disk model having magnetic field azimuthally parallel to the disk, and 

^. \ magnetic field perpendicular to the disk. Outflows with velocities, at least, ~10 

kms" 1 are driven by the protoplanets in both magnetized disk models, while 

p-) ', outflow does not appear in unmagnetized disk model. Tube-like outflows along 

the azimuthal direction of the protoplanetary disk appear in the model with 

Q^ \ magnetic field being parallel to the disk. In this model, the magnetically dom- 

I£2 [ mated regions (i.e., density gap) are clearly contrasted from other regions and 

spiral waves appear near the protoplanet. On the other hand, in the model with 
magnetic field being perpendicular to the disk, outflows are driven by a proto- 
planet with cone-like structure just as seen in the outflow driven by a protostar. 
Magnetic field lines are strongly twisted near the protoplanet and the outflows 
have well-collimated structures in this model. These outflows can be landmarks 
for searching exo-protoplanets in their formation stages. Our results indicate 
that the accretion rate onto the protoplanet tend to have a larger value than 

. !_, \ that expected from previous hydrodynamical calculations, since a fraction of the 

angular momentum of circum-planetary disk is removed by outflows, enhanced 
non-axisymmetric patterns caused by magnetic field, and magnetic braking. Pos- 
sible implications for observation are also briefly discussed. 

Subject headings: ISM: jets and outflows — MHD — stars: planetary system - 
solar system: formation 
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Introduction 



Can the giant gas protoplanets drive outflows through their formation process ? Out- 
flows are ubiquitous in star-forming region and believed to be a universal phenomenon in 
star formation process. Recently, an outflow driven by a young brown dwarf was found in 
p Ophiuchi star-forming region (Whelan et al. 2005). Recent numerical simulations show 
that protostellar outflows can be driven in the collapsing clouds when the clouds are ini- 
tially magnetized and rotating (Tomisaka 2002; Matsumoto & Tomisaka 2004; Machida et 
al. 2004, 2005b; Banerjee & Pudritz 2006). In the planet formation process, the angular 
momentum is acquired from shearing motion between the protoplanet and protoplanetary 
disk. Thus, outflows may be driven by the protoplanet, when the protoplanetary disk is 
magnetized (Sano et al. 2000; Salmeron et al. 2005; Inutsuka & Sano 2005; Matsumura & 
Pudritz 2006). Tidal interaction between the protoplanet and the protoplanetary disk tends 
to open up a gap within the protoplanetary disk where the density is so small that sufficient 
ionization is available even only by cosmic rays. The possibility to launch outflow from the 
close vicinity of a protoplanet is investigated by Fendt (2003). In his analytical study, he 
assumes the mass accretion rate and angular velocity of the protoplanet according to the 
recent results of hydrodynamical simulations, and indicates that the protoplanet can drive 
outflow with ^ 60 kms -1 of the velocity if circum-planetary disk is magnetized. However, the 
magnetic field also affects the efficiency of angular momentum transfer within and from the 
circum-planetary disk, and hence changes the mass accretion rate. Thus, three-dimensional 
MHD simulations are needed to clarify whether the proto-planet drive outflows or not. 

The giant gas planet formation has been well investigated by one-dimensional hydrody- 
namical simulations. Those studies suggest that the giant planets are formed by the gradual 
buildup of a rocky core through coagulation of solids, and this core undergoes significant gas 
accretion once the core mass reaches ~ 10 Earth mass (Mizuno 1980; Bodenheimer & Pollack 
1986; Ikoma et al. 2000). Tanigawa & Watanabe (2002a) have studied the gravitational in- 
teraction between a protoplanet and an isothermal gaseous disk using their two-dimensional 
local hydrodynamical simulation with shearing sheet model, and shown the gas flow pat- 
tern around the protoplanet. Miyoshi et al. (1999) and Tanigawa & Watanabe (2002b) 
performed three-dimensional local hydrodynamical simulation and found the detailed flow 
pattern around the protoplanet. Pioneering work on the effect of the magnetic field in proto- 
planet formation was done by Nelson & Papaloizou (2003). They have shown that accretion 
rate onto forming gas giant planet and flow pattern are modified by magnetic connection 
between the circumstellar disk and circum-planetary disk. 

In this paper, we calculate the formation of the gas giant protoplanet in the magnetized 
protoplanetary disk using three-dimensional nested grid method. We adopt the shearing 
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sheet (i.e., local) model, and found the gas giant protoplanets can drive outflows. 



2. Model and Numerical Method 

We consider only a local region around the protoplanet using shearing sheet model 
(e.g., Goldreich & Lynden-Bell 1965). We assume that the temperature is constant and 
the self-gravity of the disk is negligible in this local region. The orbit of the protoplanet is 
assumed to be circular on the equatorial plane of the disk. In our models, we set up local 
Cartesian coordinates with origin at the protoplanet's position and the x-, y-, and z-axis 
are, respectively, radial, azimuthal, and vertical direction of the disk. We solve ideal MHD 
equation without self-gravity: 

^ + V-( P v) = 0, (1) 

dv I 

+ (t,.V)t> = — VP-Viff-20 p (lxt)), (2) 

— = V x (v x B), (3) 

where p, v, P, ip e g, Q p , z, and B are the gas density, velocity, gas pressure, effective potential, 
Keplerian angular velocity of the protoplanet, a unit vector directed to the z-axis, and the 
magnetic flux density, respectively. In the above equations, the curvature terms are neglected. 
We adopt an isothermal equation of state, P = c^p, where c s is sound speed. Because we 
neglect the self-gravity of the disk in equation (2), density is scalable in basic time evolution 
equations (1),(2), and (3). Therefore we present our result for density distribution in units 
of the initial density at the mid-plane of the protoplanetary disk. The Keplerian angular 
velocity of the protoplanet is given by Cl p = (GMc/rp 1 ) 1 / 2 , where G, M c , and r p are the 
gravitational constant, mass of the central star, and the distance between the protoplanet 
and the central star, respectively. Our calculations are normalized by unit time, fip 1 , unit 
velocity, c s , and unit length, h = c s /Qp. The unit length corresponds to the scale height 
of the disk. The effective potential ip c g in our normalization is given by ip e g = — (3x 2 — 
z 2 )/2 — 3rn/(r 2 + e 2 ) 1 / 2 . The first term is composed of the gravitational potential of the 
central star and the centrifugal potential, and higher orders in x, y and z are neglected. 
The second term is the gravitational potential of the protoplanet, where rn, r, and e are 
Hill radius, the distance from the center of the protoplanet and softening parameter. The 
Hill radius are defined by rn = (M p /3M c ) 1 / 3 r p , where M p is the mass of the protoplanet. 
In the unmagnetized disk, gas flow is characterized by only one parameter, rn- In this 
paper, we adopt rn = h. This corresponds to the actual mass of protoplanet related to the 
orbital radius and the mass of the central star as M p = 3M c (h/r p ) 3 M & . In the standard 
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solar nebular (Hayashi et al. 1985), the scale height of the disk at r p = 5.2 AU is h — 
0.262 AU. Thus, the mass of the protoplanet is M p = OAMj for M c = 1M , where Mj 
is the mass of Jupiter. We will discuss our results assuming r p = 5.2 AU and M c = 1M Q 
in the following. We will show the results with other sets of parameters in a subsequent 
paper. The gas flow has a constant shear in the x-direction as Vq = (0, —3/2 a;, 0). The 
density is given by p = a exp(—z 2 /2h 2 )/(27rh 2 ) 1 ^ 2 , where a (= j°° pdz) is the surface 
density of the unperturbed disk. The density is normalized by po — &o/h (for detail see 
Miyoshi et al. 1999). We calculate three models: one unmagnetized disk model (Model 
BO), and two magnetized disk models (Models BZ and BY). In model BZ, the magnetic 
field is perpendicular to the disk with uniform field strength as B = (0, 0, 0.1 B c ), where we 
adopt central magnetic field as B c = {%ixc 2 po) l l 2 . On the other hand, the magnetic field in 
Model BY is parallel to the azimuthal direction of the disk and decrease with increasing z as 
B = (0, B c exp[—z 2 /2h 2 }, 0). In this paper, we show only strongly magnetized disk models 
in which the magnetic field energies are comparable to the thermal one, in order to see the 
magnetic effect. 

We adopt the nested grid method (for detail, see Machida et al. 2005a) to obtain high 
spatial resolution near the origin. Each level of rectangular grid has the same number of cells 
(= 64 x 128 x 32), but cell width As(l) depends on the grid level /. The cell width is reduced 
1/2 with increasing grid level (I —>■ I + 1). We use 6 grid levels (I =1,2- • -6). The box size 
of the coarsest grid I — 1 is chosen to (x, y, z) = (16h, 32h, 8h), and that of the finest grid is 
(x,y,z) = (0.5h,h,0.25h). The cell width of the coarsest grid is As(l) = 0.25/i, while that 
of the finest grid has As(6) = 7.8 x 10~ 3 h. We assume the fixed boundary condition in the 
x- and z-direction and periodic boundary condition in the y-direction. We adopt softening 
length as e = 0.03h to mimic the object with a finite size. 



3. Results 

Figure 1 shows the density distributions and velocity vectors at t — 6.11 for model B0. 
We overplot 5 or 3 levels of grids. The protoplanet revolves counterclockwise direction around 
the central star approximately once by this time. In our calculations, the flow patterns in all 
models seem to become steady on the equatorial plane in t ~ 5. Although we calculate the 
evolution until t ~ 10, the flow patterns hardly change during 5 <^t <^ 10. In Figure 1 a, the 
gas flows according to the Keplerian shear motion outside the Hill radius (r > 1): gas enters 
from upper y boundary for x > (from lower y boundary for x < 0) and goes downward 
(upward for x < 0). In the region of \x\ ^ rn, the flow is directed to the protoplanet, and 
then the gas is confronted by a shock wave as shown in Figure 1 a. Figure 1 b is enlargement 
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of Figure 1 a. In this panel, the gas flow is bent due to shock wave. Then, gas turns round by 
the Coriolis force and goes upward in the region of x < 0. When the post-shocked gas flows 
into r <^ 0.5, the gas is accreted to the circumplanetary disk. The disk rotates circularly 
around protoplanet in the counterclockwise direction near the protoplanet (r <^ 0.5). In 
Figure 1 c, the shock wave is observed on the vertical plane. Once the gas is vertically 
rolling up, it falls to the center inside the shock region. 

Figure 2 shows the density and velocity distributions for model BZ. Although the large 
scale structure in the z = plane for model BZ (Fig. 2 a) are similar to that for model BO 
(Fig. 1 a), ellipsoidal density gaps appear at (x, y) ~ (±1-5, =f3) that correspond to the 
low beta (magnetically dominated) regions. In Figure 2 a and 26, the regions enclosed by 
the red contour lines indicate magnetically dominated region (j3 p < 1), where (3 P denotes 
f3 p = 8irc^p/B 2 . In this region, the magnetic pressure exceeds the thermal pressure, while 
the thermal pressure dominates outside this region. The low density regions appear also 
inside the Hill radius (r < 1) that also correspond to the low beta regions (Fig. 2 b). In the 
neighborhood of protoplanet (r <^ 0.5), the gas enters from lower left and upper right and 
flow out to lower right and upper left, as seen in Figure 2 b. Thus, flow patterns in model 
BZ are considerably different from those in model BO near the protoplanet. The structure of 
vertical direction in model BZ is also considerably different from that in model BO. Outflow 
has appeared in model BZ since t ~ 0.66, while outflow has not appeared in model B0. White 
thick line in Figure 2 c and d mean the boundary between the outflow and inflow (contour 
curve of v z =0). We can see from Figure 2 c and d that outflow speed is much faster than the 
infall speed. Therefore, gas density decrease steeply along the z-axis as shown in Figure 2 c 
and d, because strong outflow vertically blows off the gas. The outflow took the maximum 
speed of v out = 14.1 at t — 3.37. Then, outflow became roughly steady state with velocity of 
f out ~ 10 and gradually increases again in stages near the end of the calculation. Magnetic 
field lines are strongly twisted near the protoplanet, and outflow has a well-collimated (or 
cone-like) structure that is similar to the outflow driven by a protostar (Tomisaka 2002). 

An outflow also appears in model BY. Cut planes for model BY are shown in Figure 3. 
Comparing Figure 3 a with Figure 1 a, the gas flow in model BY is similar to that in 
model B0 on the large-scale. However, the flows in small-scale are different as indicated 
by Figure 36 and lb. Spiral waves are appeared in model BY (Fig. 3 b), while the central 
region has a round shape in model B0 (Fig. 1 b). These patterns seem to be caused by 
magnetic force. The magnetic and thermal pressure dominated regions are clearly separated 
from each other in this model. The gaps in Figure 3 a and b correspond to the magnetic 
pressure dominated region as well as those of model BZ. The gas flows out of this region 
because of the strong magnetic pressure gradient. Outflow has appeared since t ~ 3.6 in 
model BY. Inside the white thick lines that correspond to the outflow region, we can see the 
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density contours expand outside, which indicates the mass ejected from the vicinity of the 
protoplanet owing to the outflow. The maximum speed of the outflow changes with time. 
The outflow becomes weak at t — 4.65 and has a i> out = 6.03 of the maximum outflow speed 
at this epoch. Then outflow becomes strong again, and has v max = 12.8 at the end of the 
calculation (t = 13.4). We found that outflows have appeared just around low beta regions 
and the outflowing region has a hollow tube structure as discussed in Fendt (2003). 



4. Discussion 

As shown in Figure 2 and 3 the gas is flowing away from the Hill region despite the 
gravity from the protoplanet. The escape velocity at r = rn = 0.04r p is 1.61 km s" 1 which 
is slower than the escape velocity from the central star 18.5 km s^ 1 in our settings. The 
sound speed at the orbit of Jupiter (5.2 AU) corresponds to c s = 0.66kms _1 in the standard 
solar nebula (Hayashi et al. 1985). Thus the outflow speeds at the end of the calculations 
correspond to 9.31kms _1 for model BZ and 8.45 kms -1 for model BY, if the protoplanet 
resides at 5.2AU in our models. The escape velocity from the protoplanet (1.61 km s^ 1 ) at the 
Hill radius is slower than the outflow speed derived from the final snapshot of our calculation. 
As shown in Figure 2 and 3, the outflow region expands far beyond the Hill radius (r ^> 1). 
In a region far from the protoplanet, outflow has the velocity much faster than the escape 
velocity from the protoplanet. If the outflow speed from the protoplanet exceeds the escape 
velocity of the central protostar, gas is continuously accelerated and eventually ejected into 
the ambient interstellar medium. The escape velocity of the central star (18.5 km s -1 ) at the 
orbit of Jupiter is about two times faster than outflow speed derived from the final snapshot 
of our calculations. However, the outflow speed for model BZ continues to increase until the 
end of our present calculation. Outflows may increase their speed with increasing mass of 
the protoplanet, while the mass of the protoplanet is assumed to be fixed in our calculations 
for simplicity. If the outflow driven by the protoplanet has the speed larger than the escape 
velocity of the central star, we can observe it analogously with the protostellar outflow. Thus 
we will possibly infer the existence of the protoplanet from the observed outflow driven by 
the protoplanet with future high resolution observational facility such as ALMA. 

Outflow may influence the evolution of the protoplanet. Angular momentum is signifi- 
cantly transported by the outflow as in the star-formation process (Tomisaka 2002). Angular 
momentum is also transported by non-axisymmetric pattern near the protoplanet. In our 
calculations, the spiral patterns are formed near the protoplanet (r <^ 0.5) in magnetized 
disk models as seen in Figure 2 b and Figure 3 b, while no non-axisymmetric patterns ap- 
pears within r <^ 0.5 in unmagnetized disk model (Fig. 1 b). Furthermore, the effect of 
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magnetic braking is important for transportation of the angular momentum. To illustrate 
this, we measured the mass and angular momentum of the gas inside r < rjj at the end of 
the calculation, and found that the acquired angular momenta per accreted mass are 0.083 
(model BO), 0.058 (model BY), and 0.062 (model BZ) in units of (c§/ft p ). Note that we 
followed only later phase of the giant planet formation; the initial mass of the protoplanet is 
already 0.4 Mj. Nonetheless, our result shows that the angular momentum of the gas accret- 
ing onto the protoplanet is more effectively transferred outside in magnetized disk models. 
The evolution of the protoplanet's rotation should be studied in long-term calculations. 
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Fig. 1. — The density (color-scale and contour) and velocity distributions (arrows) on the 
cross-section in z = (equatorial plane; a and b), y = (radial direction; c), and x = 
plane (azimuthal direction; d) for model BO are plotted at t ~ 6.11. Three levels of grid 
(I = 4, 5, 6) are overplotted in panel b, while five levels of grid (I = 2, • • •, 6) are oveplotted 
in panels a, c and d. The scale of the velocity vectors, and grid level (/) are shown in each 
panel. The unit of the velocity (i.e., sound speed) corresponds to 0.66 kms -1 , when we 
assume the orbit radius, gas temperature and mean molecular weight as 5.2AU, 122.7 K, 
and 2.34, respectively. 
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Fig. 2. — Same as Figure 1 but for model BZ at t cz 6.05. Thick red line in panels a and b 
means the magnetically dominated region (/3 P < 1, where j3 p is plasma beta). Thick white 
lines in panels c and d are the boundary between the inflow and outflow region (contour line 
of v z = 0). Inside white line, gas is outflowing. 
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Fig. 3. — Same as Figure 2 but for model BY at t ~ 6.1. 



